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摘要 :以 流 河 流域 为 研究 区 域 ,基于 SW AT 模型 模拟 1976—2012 年 洪 河 流域 的 水 沙 过 程 ,分 析 2000 年 京 津 风沙 源 治理 项 目 实 
施 前 后 流域 产 水 产 沙 时 空格 局 变化 ;研究 1980—2010 年 流域 景观 格局 变化 特征 ,揭示 景观 格局 变化 的 水 沙 响应 ;应 用 Spearman 
相关 分 析 法 分 析 流 域 景观 格局 变化 对 水 沙 过 程 的 影响 。 结 果 表 明 :与 20 世纪 80 年 代 相 比 ,2010 年 流域 林地 和 建设 用 地 增加 ， 
其 他 用 地 类 型 减 小 ;流域 景观 集中 程度 提高 .连通 性 变 优优 势 斑 块 显 著 .形状 趋 于 规则 多样 性 减少 破碎 化 程度 降低 .景观 类 
型 向 非 均衡 方向 发 展 ;流域 年 均 地 表 径 流 减 少 9mm, 产 水 量 增 加 5.44mm, 产 沙 减 小 1.59vhm”; 地表 径流 减少 区 域 占 全 流域 
89.32% , 产 水 量 增 加 区 域 占 76.71% , 产 沙 量 减 少 区 域 占 93.89% ;地 表 径 流产 水 . 产 沙 与 林地 面积 呈 负 相关 , 产 水 与 草地 面积 呈 
正 相关 ,地 表 径 流 . 产 沙 与 农业 用 地 面积 呈正 相关 ;地 表 径 流 . 产 水 . 产 沙 与 景观 形状 .Shannon's 均匀 度 .景观 分 离 度 呈 正 相 关 ， 
与 草 延 度 .最 大 斑 块 指数 呈 负 相关 ; 产 水 . 产 沙 与 斑 块 密度 和 Shannon's 多 样 性 指数 呈正 相关 ;工程 治理 后 ,流域 年 均 径 流量 与 产 
沙 量 显著 下 降 , 产 水 产 沙 高 值 区 显著 缩小 , 产 沙 关键 区 域 仍 需 治 理 。 
关键 词 :景观 格局 ; 产 沙 ;SWAT 模型 ;时 空 分 布 ; 深 河 流域 
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Abstract: The Luanhe River Basin is one of the most important sources of sand for sandstorms in Beijing and Tianjin. To 
alleviate the hazards, a sandification control program for areas near Beijing and Tianjin was carried out in June 2000, which 
could have changed the landscape pattern of this basin. Firstly, the present study used landscape pattern indices to analyze 
landscape pattern changes from 1980 to 2010. Secondly, the SWAT model was used to reveal how the changes in landscape 
pattern influenced the surface runoff, water yield, and sediment yield in 1980 and 2010, and to evaluate the effects of the 
sandification control program. Additionally, the Spearman Correlation was used to analyze the relationship between 
landscape pattern changes and flow-sediment process changes. Furthermore, the SWAT model was applied to study the 
spatial and temporal distribution of the water and sediment yield in the basin before the implementation of the program 
(1976—1999) and after the implementation of the program ( 2000— 2012) , respectively, and identify the critical zones of 


sediment yield that should be controlled at present in the Luanhe River Basin. The results were as follows; (1) compared 
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with the area of landscape types in 1980s, forest and construction land increased obviously in 2010, whereas the other land 
use types decreased in 2010. Compared with the landscape pattern in 1980, increases of splitting index ( SPLIT), 
Contagion Index ( CONTAG) , Landscape Shape Index ( LSI) , and Largest Patch Index ( LPI) , and decreases of Shannon’ 
s Diversity Index ( SHDI) and Shannon's Evenness Index ( SHEI) , diversity of landscape patterns and fragmentation, and 
more regular landscape shapes and imbalanced development of landscape types in 2010 were also observed. (2) Compared 
with the flow and sediment processes in 1980, annual average surface runoff was significantly reduced in 2010, with a 9 mm 
decrease, and the decrease in surface runoff accounted for 89.3296 of the whole basin; the annual average water yield 
increased strongly in 2010, with a 5.44mm increase, and areas of the increase in water yield accounted for 76.7196 of the 
basin; average annual sediment yield was obviously reduced, with a decrease of 1.59 t/hm^ , and 93.8996 of the basin area 
showed this trend. (3) Variations in landscape indices were correlated with the processes of flow and sediment well. The 
surface runoff, water yield, and sediment yield were negatively correlated with the area of forest. The water yield was 
positively correlated with the area of grassland. The surface runoff and sediment yield were positively dependent on the area 
of agricultural land. Surface runoff, water yield, and sediment yield were positively correlated with SPLIT, LSI, and SHEI. 
On the contrary, surface runoff, water yield, and sediment yield were negatively correlated with LPI and CONTAC. 
Additionally, water and sediment yield both positively correlated with SHDI ( Shannon's Diversity Index) and PD ( Patch 
Density). (4) Soil erosion in most zones of Luanhe River Basin has been reduced with the implementation of the program. 


However, a few zones were still experiencing severe soil erosion that should be taken into account. 


Key Words: landscape pattern; sediment yield; SWAT model; temporal-spatial distribution; Luanhe River Basin 


景观 格局 是 不 同 景 观 馈 黄体 在 空间 上 的 排列 ,是 不 同 生 态 系 统 组 成 的 地 理 综合 体 ,能 有 效 揭示 区 域 生 
态 状 况 和 空间 变异 特征 :2 。 景 观 格局 变化 主要 是 人 类 活动 通过 土地 利用 变化 驱动 的 ?1 。 景 观 的 空间 配置 与 
类 型 组 成 影响 水 循环 ,从 而 影响 流域 系统 的 水 文 过 程 “; 。 从 景观 生态 学 角度 探讨 流域 景观 变化 的 水 沙 过 程 
响应 ,能 够 揭示 人 类 活动 对 流域 水 沙 的 影响 。 在 水 资源 短缺 和 水 土 流失 严重 的 区 域 , 开 展 此 研究 ,可 以 为 区 域 
土地 利用 的 合理 调整 .水土 保持 治理 与 生态 恢复 提供 科学 依据 。 目 前 ,景观 格局 变化 的 水 沙 过 程 响应 已 成 为 
研究 热点 。Dixon 5j Earls! 发 现 佛罗里达 州 查理 溪流 域 径流 随 城市 扩张 而 增加 ; 王 计 平等 发 现 斑 块 密度 是 影 
响 黄土 丘陵 沟 蜜 区 产 沙 的 重要 指标 1 ; 黄 志 霖 等 发 现 三 峡 库 区 景观 格局 对 径流 和 泥 沙 影响 均 显 著 ! 7 ; 郭 军 庭 
等 发 现 潮 河 流域 灌木 林地 和 草地 能 够 增加 径流 ; 林 炳 青 等 发 现 晋江 流域 洪水 径流 与 景观 格局 相关 性 显 
AU ; 李 唱 与 周 自 翔 发 现 延 河 流域 景观 格局 变化 与 产 沙 显著 相关 "” ; 郝 振 纯 等 分 析 土 地 利用 情景 变化 对 海河 
流域 典型 区 域 径流 的 影响 。 

流 河 流域 距 京 津 地 区 约 240km ,是 京 津 沙尘暴 的 重要 通道 和 沙 源 地 之 一 号 。 北 京 沙 尘 绝 大 部 分 源 自 外 
地 ,与 洪 河 流域 的 丰 宁 和 多 伦 的 土壤 属性 非常 相似 。 为 优化 京 津 地 区 环境 质量 ,2000 年 京 津 风沙 源 治理 
工程 启动 ,多 伦 县 为 工程 发 端 地 !21 。 淋 河流 域 位 于 工程 治理 区 内 :5 ,主要 通过 土地 利用 变化 影响 地 表 水 下 
18 . 蒸 散 (发 ) 等 水 文 过 程 , 因 此 ,探讨 景观 格局 变化 的 水 文 效 应 ,对 流 河 流域 水 土 流失 治理 具有 重要 意义 。 目 
前 , 流 河 流域 产 流产 沙 相关 研究 得 到 广泛 关注 ,研究 多 集中 于 径流 ,研究 范围 主要 为 潘 家 口水 库 上 游 :5 与 全 
流域 59 ; 泥 沙 研究 主要 针对 潘 家 口水 库 上 游 流域 ' 引 , 较 少 涉及 全 流域 ;水 文 影响 因素 主要 集中 于 土地 利用 与 
气候 变化 对 径流 的 影响 , 较 少 从 景观 格局 变化 分 析 产 流产 沙 。 

本 文 运用 景观 格局 指数 分 析 1980—2010 年 洪 河流 域 景观 格局 变化 特征 ,结合 SW AT 模型 探讨 景观 格局 
变化 的 流域 水 沙 过 程 响应 ;应 用 Spearman 相关 分 析 法 分 析 景 观 格局 变化 对 水 沙 过 程 的 影响 ;应 用 SW AT 模型 
模拟 1976—2012 年 深 河 流域 径流 与 泥 沙 ,分 析 流 域 2000 年 治理 工程 实施 前 后 多 年 年 均 产 流 与 产 沙 的 时 空格 
局 变化 ,识别 流域 仍 需 重点 治理 的 产 沙 关键 源 区 ,以 期 为 洪 河 流域 的 后 续 治理 与 土地 利用 规划 提供 科学 依据 。 
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1 研究 区 域 ,方法 与 数据 来 源 


1.1 研究 区 域 

漆 河 流域 位 于 华北 北部 , 深 河 全 长 888km ,流域 总 面积 为 44750km”( 图 1)。 深 河 发 源 于 河北 丰 宁 , 流 经 内 
蒙古 ,河北 及 辽宁 的 7 市 ( 盟 ) 的 27 个 县 (区 )' 沾 。 漆 河流 域 属于 典型 温带 大 陆 性 季风 气候 ,年 均 气 温 5 一 
12%C ,多 年 平均 降水 量 为 400— 700mm ,主要 集中 分 布 于 7 一 8 月 。 漆 河流 域 地 形 起 伏 较 大 ,由 西北 向 东南 倾 
斜 ,分 为 坝 上 高 原 .燕山 山地 南部 平原 3 种 地 貌 类 型 。 流 域内 主要 土地 利用 类 型 为 林地 、 草 地 和 农业 用 地 。 
流域 土壤 类 型 丰富 ,主要 包括 草 和 旬 土 .森林 土 REE RR 潮 士 和 褐 土 等 '” 。 
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图 1 漆 河 流域 位 置 图 
Fig.1 Location of Luanhe River Basin 


1.2 ”数据 来 源 与 处 理 

DEM :90m 分 辨 率 , 来 源 于 美国 地 质 勘 探 局 (USGS) 。 根 据 流域 情况 与 相关 研究 , 设 定 最 小 集 水 面积 
IB] (ÉL 20000ha ,流域 划分 为 116 个 子 流 域 ,以 子 流域 为 基础 单元 研究 产 水 产 沙 的 空间 分 布 特征 。 根 据 土 壤 、 
土地 利用 坡度 的 分 布 特征 ,面积 阔 值 分 别 设 定 为 13% 10% .15% ,生成 1521 个 HRU ,充分 体现 了 洪 河 流域 下 
热 面 的 空间 不 均匀 性 。 

土地 利用 数据 :1980s 土地 利用 图 ,分辨 率 1km ,来 源 于 地 球 系统 科学 数据 共享 平台 ;2010 年 土地 利用 数 
据 源 于 分 辩 率 30m 的 TM 数据 (轨道 号 为 121/32,121/33,122/31,122/32,123/31,124/31) ,来 源 于 USGS., 
运用 ENVI4.8 和 ArcGIS 9.3 £e Bé .拼接 等 处 理 , 采 用 人 机 交互 解 译 影像 ,并重 采 样 为 Ikm 分 辨 率 。 依 据 
土地 资源 分 类 系统 ,采用 国家 标准 体系 ,分 为 林地 、 草 地 农业 用 地 ,水域 .建设 用 地 和 未 利用 地 六 类 。 依 据 
数据 可 获得 性 研究 目的 以 及 20 世纪 80 到 90 年 代 土 地 利用 变化 相对 缓慢 '2 ,采用 1980s 土地 利用 数据 用 于 
1976—1999 年 时 间 段 的 模拟 ; 同样 地 ,2010 年 土地 利用 数据 用 于 2000—2012 年 时 间 段 的 模拟 。 

土壤 数据 :1:100 万 的 世界 土壤 数据 库 (HWSD ) 来 源 于 寒 区 旱 区 科学 数据 中 心 。 土 壤 水 文 分 组 根据 在 完 
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全 湿润 并 不 冻 的 条 件 下 土壤 最 小 渗透 属性 来 确定 '” ;土壤 容重 等 采用 SPAW 软件 进行 计算 。 

气象 数据 :来 源 于 中 国 气象 科学 数据 共享 服务 网 ,包括 多 伦 、 丰 宁 、 围 场 、 承 德 \ 乐 亭 5 个 气象 站 1961 一 
2012 年 的 日 降水 气温、 风速 等 。 利 用 上 述 气 象 数据 建 立 天 气 发 生 器 ,输入 数据 由 pcpSWAT 和 dewSWAT 软 
件 得 出 5 。 

水 文 数据 :1976 一 1991 年 承德 .下 板 城 . 李 营 、 宽 城 . 漆 县 5 个 水 文 站 的 月 径流 量 潘 家 口 和 大 黑 洒 水库 月 
出 流 数据 和 漆 县 站 泥 沙 数据 ,来 源 于 深 河 流域 水 文 年 鉴 等 资料 。 

SWAT 模型 空间 数据 统一 采用 Albers 等 积 投 影 。 借 鉴 相关 研究 "1 ,90m 分 辩 率 DEM 满足 流域 水 系 的 
正确 划分 ;1km 分 辨 率 土 地 利用 图 与 土壤 图 适用 于 流域 水 沙 过 程 模拟 。 

1.8 ”研究 方法 
1.3.1 SWAT 模型 

SWAT 模型 的 开发 目的 是 模拟 与 预测 具有 复杂 变化 的 土壤 、 土 地 利用 和 管理 方式 的 较 大 流域 中 土地 利用 
及 管理 方式 对 流域 产 流 、 产 沙 及 污染 物 负 荷 的 影响 。 文 中 选用 SCS 径流 曲线 数 方 法 计算 地 表 径 流 , 泥 沙 由 改 
itf] MUSLE 方法 计算 ,采用 Penman-Monteirh 方法 计算 潜在 蒸发 量 , 马 斯 京 根 方法 进行 河道 演算 。 采 用 
SWAT2009 的 敏感 性 分 析 模 块 中 全 局 参数 敏感 性 分 析 方 法 (LH-OAT) 进行 径流 和 泥 沙 参数 敏感 性 分 析 '*] 。 
运用 SWAT-CUP 软件 的 t 检 验方 法 进行 参数 的 深层 筛选 六 ;依次 校准 径流 与 泥 沙 ,通过 SWAT-CUP 的 SUIF- 
2 方法 (站 对 径流 进行 初步 率 定 后 手动 调 参 ; 泥 沙 校准 采用 手动 调 参 。 

Nash-Suttcliffe 系数 (Ens) ,决定 系数 (R? ) 和 相对 误差 (Re) 为 模型 适用 性 评价 指标 。 当 Ens>0.5,R >0.5， 
Re<20% 时 ,径流 模拟 结果 是 满意 的 ;Re<55% 时 , 泥 沙 模拟 结果 是 满意 的 '*1。 

1974 一 1975 年 为 预 热 期 ,1976 一 1983 年 作为 模型 校准 期 ,1984 一 1991 年 为 验证 期 。 利 用 1992 一 2012 年 
气象 数据 、1980s 和 2010 年 土地 利用 数据 和 已 率 定 参 数 , 得 到 1992—2012 年 径流 与 泥 沙 模拟 数据 。 借 鉴 相 关 
WRAN :保持 气象 .土壤 等 输入 条 件 不 变 , 仅 改 变 土地 利用 条 件 时 ,水 文 模型 输出 结果 差异 仅 是 由 土地 利 
用 变化 造成 的 。 本 文 以 相同 的 气象 .地 形 .土壤 .水 库 出 流 等 条 件 ,分 别 输入 1980s 和 2010 年 土地 利用 图 , 模 
拟 和 分 析 景 观 格局 变化 的 水 沙 响应 ,评价 工程 治理 前 后 流域 水 土 流失 变化 。 

针对 上 述 假设 ,采用 2000—2012 年 气象 数据 , 满足 假设 中 相同 气象 条 件 。 研 究 区 面积 较 大 ,2000 一 2012 
年 全 流域 尺度 地 形 地 貌 变 化 土壤 类 型 变化 皆 可 忽略 ,满足 假设 中 相同 的 地 貌 与 土壤 条 件 。 漆 河 干流 建设 了 
潘 家 口 和 大 黑 汀 水库, 支流 的 小 型 水 库 对 流域 整体 水 文 过 程 影响 较 小 ,因此 主要 考虑 语 ,大 水 库 ( 分 别 由 1981 
年 .1983 年 开始 运行 至 今 "“29 ) ,满足 相同 的 水 利 工 程 条 件 。 可 见 假 设 合理 ,研究 方法 可 行 。 

1.3.2 ”景观 格局 分 析 

运用 ENVI 4.8 计算 景观 类 型 面积 变化 特征 ;利用 FRAGSTATS 3.3 软件 分 析 流 域 景观 格局 指数 变化 ;应 用 
Spearman 相关 分 析 法 六 ,分 析 流 域 水 沙 过 程 变化 与 景观 格局 指数 变化 的 相关 关系 。 

考虑 相关 文献 中 对 水 沙 影响 较 大 、 意 义 明确 且 相 对 独立 的 景观 指数 ' ,本 文选 取 斑 块 类 型 面积 (CA) JE 
块 密 度 (PD)、 景 观 形状 指数 (LSI)、Shannon 多 样 性 指数 (SHDI)、Shannon 均匀 度 指数 (SHEI)、 分 离 度 
(SPLIT) .最 大 斑 块 指数 (LPI) 和 草 延 度 (CONTAG ) ,分 别 反 映 景观 类 型 面积 .形状 复杂 度 、 破 碎 化 程度 、 多 样 
性 与 连通 性 。 

统计 1980 年 .2010 年 全 流域 和 各 子 流域 景观 指数 ,并 计算 全 流域 和 各 子 流域 1980—2010 年 景观 格局 的 
变化 率 , 计 算 公 式 为 : 


F. - F. 
P, = i,2010 i585. 100g, (1) 
F, 1080 
式 中 BO i KIAK TREER, F, 2010 为 2010 年 景观 格局 下 i 类 景观 格局 指数 ,FF 1 1980s 景观 格局 


下 i 类 景观 格局 指数 。 
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2 结果 与 分 析 


2.1 模型 的 校准 与 验证 

如 表 1 所 示 , 各 水 文 站 控制 的 子 流域 中 敏感 性 因子 对 径流 的 敏感 度 存在 差异 ,这 与 研究 区 面积 较 大 ,地 
E .土地 利用 与 土壤 类 型 多 样 有 关 。 然 而 在 所 有 子 流 域 中 ,径流 曲线 系数 CN2 的 敏感 度 均 较 高 ,这 与 其 他 研 
究 相 一 致 “1 ,此 参数 直接 影响 径流 的 大 小 , 值 越 高 径流 量 越 大 。 

模型 校准 与 验证 结果 见 图 2、 表 2,Ens 和 尼 均 达到 0.7 以 上 。 径流 模拟 效果 较 好 , 泥 沙 模拟 与 实测 值 的 

过 程 线 也 趋 于 一 致 ,但 丰 水 年 的 泥 沙 a 村 别 是 洪水 期 。 这 主要 由 于 SWAT 模型 在 产 沙 量 高 
于 2000t 时 会 产生 低 于 实际 值 的 模拟 结果 '”| ,而 深 河流 域 年 均 产 沙 量 i 508 00 E EE 
水 时 段 模拟 值 也 会 低 于 实测 值 '3 ea 

综 上 ,SWAT 模型 模拟 径流 与 泥 沙 的 精度 均 达 到 要 求 。 流 河流 域 面 积 较 大 ,研究 区 上 游 的 漆 河 源头 为 冬 
春季 的 风蚀 与 夏季 的 水 蚀 互 为 动力 ,但 水 蚀 仍 是 全 流域 最 主要 的 土壤 侵 刨 过程,SWAT 模型 适用 于 土壤 、 
土地 利用 等 具有 复杂 变化 的 漆 河 流域 水 沙 模拟 研究 。 
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图 2 承德 站 、 下 板 城 站 、 李 营 站 、 宽 城 站 、 漆 县 站 月 径流 与 深 县 站 (f) 月 泥 沙 模拟 
Fig.2 Monthly simulated and observed runoff for: the Chengde station, the Xiabancheng station, the Liying station, the Kuancheng 


station, the Luanxian station; monthly simulated and observed sediment for the Luanxian station 
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表 2 径流 和 泥 沙 的 校准 与 验证 结果 


Table 2 Calibrated and validated results of runoff and sediment 


类 型 水 文 站 校准 期 Calibration period 验证 期 Validation period 
Type Hydrologic station Ens R? Re /% Ens R? Re /% 
径流 Runoff 承德 0.73 0.73 8.20 0.72 0.74 -3.45 
下 板 城 0.87 0.86 1.19 0.73 0.73 -0.72 
李 营 0.81 0.82 12.27 0.80 0.87 0.80 
宽 城 0.95 0.95 -10.01 0.75 0.85 -19.38 
WH 0.94 0.92 5.22 0.74 0.80 14.09 
泥 沙 Sediment UB 0.74 0.93 40.41 0.75 0.71 -40.32 


T 


2.20 治理 工程 前 后 产 流 与 产 沙 的 时 空 分 布 特征 
2.2.1 时 间 分 布 特征 

2000 年 京 津 风 沙 源 治 理 以 来 ,年 均 径 流量 和 产 沙 量 比 1976—1999 年 分 别 减 少 74.26% 与 90.25% ,而 降水 
量 仅 减少 3.84% (图 3) ,表明 人 类 活动 (流域 治理 ) 是 影响 深 河 流域 径流 与 泥 沙 的 主要 因素 。 这 是 由 于 影响 径 
流量 变化 的 主要 因素 为 气候 变化 和 人 类 活动 ,气候 因素 主要 表现 为 降水 变化 ,人 类 活动 主要 为 消耗 用 水 和 土 
地 利用 变化 站。 相关 研究 也 认为 流域 治理 是 径流 量 减 少 的 重要 原因 "9 ,其 中 以 水 保 措施 改变 下 垫 面 的 影响 
量 约 为 消耗 用 水 的 2 倍 1。 

工程 开展 前 后 ,年 径流 与 泥 沙 均 主 要 集中 于 漆 河 流域 的 汛期 (6 一 9 月 ) ,7 一 8 月 径流 与 泥 沙 量 最 大 , 约 为 
年 径流 与 泥 沙 量 的 44% 和 84% ,治理 后 7 一 8 月 的 洪峰 径流 显著 下 降 , 工 程 实施 对 年 内 洪峰 流量 的 调 蓄 效果 
显著 (图 4)。 
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3 1976—2012 年 的 径流 与 泥 沙 年 际 变化 趋势 图 4 月 均 径流 与 泥 沙 的 变化 趋势 

Fig.3 Inter-annual trend of precipitation, runoff and sediment Fig.4 Distribution of mean monthly runoff and sediment in the 
from 1976 to 2012 in the Luanhe River Basin Luanhe River Basin 


2.2.02 ”空间 分 布 特征 

工程 实施 后 , 产 水 与 产 沙 的 高 值 区 显著 缩小 (图 5 .图 6) 1976—1999 年 ,下 游 多 数 区 域 的 径流 深 大 于 
150mm;2000 一 2012 年 ,只 有 位 于 下 游人 海口 附近 的 3 个 子 流域 径流 深 大 于 150mm, 而 中 上 游 出 现年 均 产 水 
量 小 于 50mm 的 区 域 ,如 洛 源 县 (34 .35 号 子 流域 ) 等 (图 5) 。 

如 图 6 所 示 ,工程 实施 前 ,流域 年 均 产 沙 模 数 大 于 10 t/hm? 的 区 域 位 于 下 游 的 青龙 县 ;工程 实施 后 ,年 均 
产 沙 模 数 大 于 10 t/hm^ 的 区 域 已 消失 ,而 2000 年 以 后 流域 中 上 游 的 大 部 分 区 域 年 均 产 沙 模 数 小 于 1 hm 。 
仍 需 重点 治理 的 区 域 主要 为 洪 河 上 游 的 洛 源 (35 号 子 流域 ) ET (33,38 号 子 流 域 ) .围场 (28 号 子 流 域 ) 与 
隆 化 (42 号 子 流域 ) ;中 游 的 漆 平 (71 号 子 流域 ) .承德 市 (71、74 号 子 流域 ) .承德 县 (75 .77 号 子 流域 ) .平泉 
(68,75 号 子 流域 ) 凌 源 (68 号 子 流域 ) ;下 游 的 迁 安 与 户 龙 (114 号 子 流域 ) 。 
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图 5 1976—1999 5 2000—2012 年 的 产 水 量 空间 分 布 
Fig.5 The spatial distribution map of simulated annual average water yield from 1976 to 1999 and during 2000 to 2012 
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6 1976—1999 ££ 5j 2000—2012 年 的 产 沙 量 空间 分 布 
Fig.6 The spatial distribution map of simulated annual average sediment yield from 1976 to 1999 and during 2000 to 2012 


2.4” 深 河流 域 治理 与 景观 格局 变化 

与 20 世纪 80 年 代 相 比 ,2010 年 林地 建设 用 地 分 别 增加 41.40% 577.82% ;农业 用 地 、 草 地 、 未 利用 地 和 
水 域 分 别 减少 26.9996 41.76% .0.50% 与 36.88%( 表 3)。 这 与 流域 治理 工程 关系 密切 ,1980 一 1989 年 ,流域 
上 游 开 展 人 工 造林 4700km? ,建设 草场 200km? 7" ;2000 年 国家 水 土 保持 重点 项 目 京 津 风 沙 源 治理 工程 启动 
后 ,采取 人 工 造 林 改良 草场 .提高 植被 覆盖 度 等 措施 建立 高 质量 的 流域 生态 防护 ,治理 成 绩 显 著 , 如 多 伦 县 造 
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PK 565km”! 中 ; 正 蓝 旗 造 林 262km?、 草 地 治理 445km? 7? ; 承德 县 造林 148 km? , 治理 252km2034 ; 丰 宁 造林 
138km? ,治理 248km*'” ;围场 退耕 还 林 665km C9 ; 深 平 县 森林 徐 盖 度 提高 15.4% , 治理 面积 950km? 等 。 

与 80 年 代 相 比 ,2010 年 流域 分 离 度 和 形状 复杂 度 显著 降低 ,景观 多 样 性 和 均匀 度 略 有 降低 ; 斑 块 密度 变 
化 不 大 ;最 大 斑 块 指数 . 划 延 度 显著 上 升 ( 表 4) ,说 明 林 草 建设 与 修复 工程 促进 流域 景观 破碎 化 降低 ,景观 分 
布 更 集中 , 斑 块 形状 更 规则 ,景观 连通 性 提高 .景观 类 型 向 非 均衡 方向 发 展 , 林 地 景观 优势 更 显著 。 


表 3 1980 一 2010 年 景观 类 型 转移 矩阵 
Table 3 Transition matrixes of different landscape types from 1980 to 2010 
农业 用 地 建设 用 地 


景观 类 型 林地 草地 水 域 Agriculture Construction 未 利用 地 2010 年 
Landscape type Forest Grassland Water land land Unused land 
林地 Forest/ km? 16006.38 6583.98 244.77 3626.65 15.54 167.22 26644.54 
草地 Grassland/km? 1237.64 4539.96 249.78 1525.47 8.13 309.87 7870.85 
水 域 Water/km? 46.32 147.65 256.08 173.34 2.03 4.29 629.71 
农业 用 地 Agriculture land/km? 1171.28 1715.23 200.15 4402.18 31.43 87.25 7607.52 
建设 用 地 Construction land/km? 189.15 220.41 17.2 635.58 120.69 22.2] 1205.3 
未 利用 地 Unused land/km? 192.26 307.75 29.72 57.16 0 205.19 792.08 
1980 年 /km? 18843.03 13514.98 997.7 10420.38 177.82 796.09 44750 
表 4 1980 一 2010 年 景观 指数 变化 分 析 
Table 4 Analysis of landscape metrics from 1980 to 2010 
景观 指数 斑 块 密度 。 最 大 斑 块 面积 。 蔓延 度 分 离 度 n 斑 抉 形状 
x Landscape metrics PD LPI CONTAG SPLIT SDT SHEI LSI 
- 1980 0.05 15.01 39.91 14.37 1.24 0.69 43.27 
2010 0.06 56.48 45.79 3.02 1.12 0.62 30.62 


PD: Patch Density; LPI; Largest Patch Index; CONTAG: Contagion Index; SPLIT; Splitting Index; SHDI: Shannon's Diversity Index; SHEI: 


Shannon's Evenness Index; LSI; Landscape Shape 


2.5 景观 格局 变化 的 水 沙 过 程 响 应 

如 图 7 所 示 , 与 20 世纪 80 年 代 相 比 ,2010 年 流域 年 均 地 表 径 流 减 少 9mm;109 个 子 流域 地 表 径 流 减 小 ， 
== 面积 约 39971km” , 占 全 流域 的 89.32%。 流域 年 均 产 水 量 增 加 5.44mm; 85 个 子 流 域 产 水 量 增加 ,面积 约 
34326.18km? , 占 全 流域 的 76.71%。 流 域 年 均 产 沙 减 小 1.59t/hm ; 104. 个 子 流 域 产 沙 减 小 , 面积 约 42015. 
37km? , 占 全 流域 的 93.89%。 可 见 , 漆 河流 域 水 土 保持 治理 能 够 很 好 地 调节 地 表 径 流 ,减少 泥 沙 。 

2.5.1 景观 格局 变化 的 地 表 径 流 响应 

与 20 世纪 80 年 代 相 比 ,2010 年 地 表 径 流 减 少 最 显著 的 区 域 主 要 位 于 上 游 16、35 号 子 流域 ,中 游 63、74 
号 子 流域 ,下 游 97 .98 101 号 等 子 流域 (图 7、 表 5)。 这 些 区域 地 表 径 流 减 少 与 林 草 植被 建设 密切 相关 。 例 如 
63.74 号 子 流域 位 于 承德 县 境内 ,实施 京 津 风沙 源 治理 工程 以 来 ,承德 县 调整 产业 结构 退耕 还 林 还 果 、 造 林 育 
林 , 全 县 森林 履 盖 率 达 到 6190 0 ,有 效 拦 著 降 水 ,减少 地 表 径 流 。 

2010 年 仅 12 个 子 流域 地 表 径 流 增 加 (图 7) 。 如 113 号 子 流 域 , 位 于 迁 安 境内 ,地表 径 流 增 加 最 显著 。 迁 
安 市 经 济 依托 于 铁 矿 资源 ,矿山 排 土 场 和 尾 矿 库 在 堆放 过 程 中 导致 绿化 率 较 低 。 虽 然 进行 综合 治理 ,但 是 由 
于 植物 种 类 与 苗 龄 选择 不 当 、 栽 植 方式 与 时 间 不 合理 ,因而 植被 覆盖 率 仍 然 较 低 ,而 且 河 道 两 侧 山体 上 的 排 土 
场 无 法 进行 清除 ,难以 采取 水 保 措 施 '3 ,水 土 流 失 傅 加 严重 ,因而 地 表 径 流 增加 。 

2.5.2 景观 格局 变化 的 产 水 响应 

与 1980 年 相 比 ,2010 年 产 水 量 增加 最 显著 区 域 主要 位 于 研究 区 上 游 2、3、6、10、13 号 子 流域 ,中 游 66、 
107,108 号 子 流域 和 下 游 106、109 .111 .113 一 115 号 子 流 域 (图 7、 表 5)。 产 水 减少 区 域 主要 位 于 深 河 上 游 的 
16,17,19 号 等 子 流域 ,中 游 的 48、51、58 号 等 子 流域 ,下 游 的 101、104、105、116 号 子 流域 。 在 产 水 量 增加 的 子 
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流域 中 ,如 19 等 子 流域 地 表 水 是 减少 的 (图 7) ,这 些 区 域 多 数 由 于 水 土 保持 治理 程度 提高 ,地 表 径 流 减少 
时 大 量 转化 补充 土壤 含水 量 与 地 下 径流 。 在 产 水 量 减少 的 子 流域 中 ,如 流域 上 游 的 围场 .河源 与 丰 宁 卉 内 
16,17,29 号 等 子 流域 地 表 水 也 是 减少 的 (图 7) ,这 是 由 于 沾 源 县 植被 覆盖 以 草地 为 主 、 丰 宁 与 围场 的 林地 建 
设 仍 需 继续 开展 ,地 表 径 流 仍 为 产 水 的 主要 组 成 部 分 , 壤 中 流 与 地 下 径流 的 增加 程度 尚 显 不 足 。 


地 表 径流 变化 /mm 
E -76—20 
ES -20—-5 
E -5 一 0 
E 0—5 
EB 5 一 14 


产 水 量变 化 /mm 产 沙 量变 化 /(t/hm?) 


EE -48—20 EB -10—5 
ES -20—-5 EH — 
E -5—0 — Mani 
ES 0 一 15 EB 0 一 1 
EB 15 一 53 mt :— 
0 50 100km 
L1. 1 


图 7 深 河 流域 1980—2010 年 年 均 地 表 径流 、 产 水 量 、 产 沙 量变 化 的 空间 分 布 
Fig.7 Spatial distribution of average yearly surface runoff, water yield and sediment yield from 1980 to 2010 


2.5.3 ”景观 格局 变化 的 产 沙 响应 

与 20 世纪 80 年 代 相 比 ,2010 年 景观 格局 下 , 产 沙 减 少 相 对 较 显 车 的 区 域 主要 位 于 漆 河 流域 上 游 的 21、 
23 .30 .35 号 子 流域 ,中 游 63 号 子 流域 和 下 游 的 105,111 号 子 流域 (图 7、 表 5)。 结 果 表 明 , 这 些 区域 水 土 保 
持 治理 效果 最 显著 ,土壤 侵蚀 程度 大 大 降低 。 例 如 63 号 子 流域 ,位 于 承德 县 东北 部 , 产 沙 量 与 地 表 径流 明显 
减少 ,这 与 承德 全 县 森林 歼 盖 率 增加 有 关 '” ,森林 覆盖 率 增加 能 够 改善 地 表 径流 .抑制 产 沙 。 
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m5 地 表 径流 、 产 水 与 产 沙 变化 的 空间 分 布 
Table 5 Spatial distribution of surface runoff, water yield and sediment yield changes 
类 型 变化 子 流域 所 在 区 域 ( 县 区、 旗 ) 
Type Chang Subbasin Location 
地 表 径 流 nm 围场 沾 源 、 承 德 市 双 桥 区 、 承 德 县 .青龙 、 建 
Surface runoff 减少 最 显著 16,35,63 74 .97.98 101 .105 112 116 B 迁 安 卢 龙 
兽 加 最 显著 79 113 承德 县 . 迁 安 
产 水 量 成 少 最 显著 35 .63 75,116 MU JRE RE JEE BR RASKE 
Water yield "ITE 2,3, 6, 10, 13, 66, 79, 106—109, 111, 正 蓝 旗 、 克 什 克 腾 旗 、 多 伦 、 承 德 市 双 桥 区 、 宽 
iia 113—115 城 迁 西 青龙 县 、 迁 安 、 卢 龙 、 漆 县 
产 沙 量 减少 最 显著 21.23.30.35.63 105 111 围场 . 沾 源 ,承德 县 青龙、 建 昌 、 迁 西 
Sediment yield 兽 加 最 显著 20.32 40.75.85 丰 宁 EHA ETE 、 宽 城 


产 沙 增加 的 子 流域 仅 12 个 ,最 显著 的 位 于 流域 上 游 的 20 .32 ,40 号 子 流域 ,中 游 75 .85 号 子 流域 (图 7、 表 


5), 2010 年 ,20 .32 40 5T Ep 


地 ,草地 减 沙 ,农业 用 地 增 沙 ,因而 水 土 流失 严重 、 产 沙 增加 。 


2.6 


景观 格局 对 水 沙 过 程 的 影响 


也 减少 被 农业 用 地 占用 ;75 、85 号 子 流 域 草 地 减少 被 农业 用 地 占用 。 林 


景观 类 型 与 水 沙 过 程 相关 性 见 表 6, 林 地 与 地 表 径 流产 水 量 和 产 沙 量 呈 负 相 关 关 系 , 说 明 深 河流 域 林 地 


面积 增加 对 径流 和 泥 沙 都 具有 积极 抑制 作用 ,这 是 由 于 森林 植被 的 下 渗 与 蒸发 


量 较 大 ,可 以 有 效 的 调节 地 表 


径流 与 洪峰 流量 ,上 且 土壤 含水 量 的 增加 也 将 减少 泥 沙 的 产生 与 输 移 '” 。 草 地 面积 与 产 水 量 呈 正 相 关 ,与 产 沙 
呈 负 相关 ,表明 草地 面积 的 增加 会 促进 产 流 ,减少 产 沙 ,主要 由 于 草地 能 将 坡 面 径流 流速 减 小 40% ,将 起 流 时 
EHER 3 倍 ,增加 径流 入 渗 率 .土壤 含水 量 和 地 下 水 ,从 而 增加 产 流 “” ;草地 的 植被 根系 存在 于 土壤 表层 , 提 
高 土壤 对 径流 侵蚀 产 沙 的 抵抗 力 。 农 业 用 地 与 地 表 径 流 . 产 沙 量 均 呈正 相关 ,由 于 农业 用 地 的 土壤 结构 被 破 
坏 , 降 低 了 径流 下 渗 速 率 .土壤 含水 量 和 地 下 水 补给 ,促使 地 表 径 流 增加 ; 而 耕地 植被 覆盖 度 较 低 ,水土 流 


失 也 相对 容易 。 建 设 用 地 与 未 利用 地 与 产 沙 相关 性 较 弱 ,可 能 与 两 者 面积 小 有 关 。 


Ro ， 湾 河流 域 景观 格局 指数 变化 与 产 流产 沙 变化 间 的 相关 性 
Table 6 The correlation coefficients among landscape metrics and water and sediment yield 
类 型 景观 格局 指数 地 表 径 流 产 水 量 产 沙 量 
Type Landscape metrics Surface runoff Water yield Sediment yield 
景观 类 型 林地 Forest 一 0.640 ** -0.505 ** -0.738 ** 
Landscape type 草地 Grassland 0.320 0.633 ** -0.452 ** 
水 域 Water 0.023 0.078 0.036 
建设 用 地 Construction land 0.029 0.119 0.203 * 
未 利用 地 Unused land 0.001 -0.005 0.184 * 
农业 用 地 Agriculture land 0.382* —0.053 0.586 ** 
景观 水 平 斑 块 密度 PD 0.184 0.323 ** 0.425 ** 
Landscape level 景观 形状 指数 LSI 0.213* 0.280 ** 0.401 ** 
最 大 斑 块 指数 LPI -0.594 * -0.627* -0.634 ** 
Shannon's 多 样 性 指标 SHDI 0.15 0.557** 0.533 ** 
Shannon's 均匀 度 指数 SHEI 0.537 * 0.524 ** 0.488 ** 
4) Ej RE SPLIT 0.414 ** 0.327 ** 0.422 ** 
草 延 度 指 数 CONTAG -0.276 ** -0.355 ** -0.254 ** 


* * 表示 相关 显著 水 习 


为 0.01( 双 尾 检测 ) , * 表示 相关 显著 水 平 为 0.05( 双 尾 检测 ) 


地 表 径 流 、 产 水 、 产 沙 与 景观 形状 复杂 程度 、Shannon's 均匀 度 、 景 观 分 离 度 呈正 相关 ， 与 蔓延 度 .最 大 斑 块 


HACE PAX ;景观 斑 块 密度 和 Shannon's 多 样 性 指数 与 产 水 产 沙 呈正 相关 ( 表 6) 。 表 明 大 力 营 造林 草地 后 ， 
景观 形状 规则 化 .连通 度 增加 流域 优势 景观 类 型 林地 增加 ,导致 流域 降水 截留 能 力 与 水 源 涵 养 能 力 增 强 ,地 
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表 径 流产 水 和 产 沙 降低 。 
3 结论 


(1) 与 80 年 代 相 比 ,2010 年 漆 河 流域 林地 、 建 设 用 地 增加 ,其 他 用 地 类 型 减 小 ;流域 斑 块 密度 变化 不 明 
显 ,景观 集中 连 片 度 提 高 .连通 性 变 优 、 优 势 斑 块 类 型 明显 .破碎 化 降低 .景观 形状 越 来 越 规则 .多样 性 程度 减 
少 ,景观 类 型 向 非 均衡 方向 发 展 。 

(2) 不 同 景观 格局 下 , 洪 河 流域 水 沙 过 程 差异 显著 。 与 1980 年 相 比 ,2010 年 流域 年 均 地 表 径 流 减 少 
9mm, 产 水 量 增加 5.44mm , 产 沙 减 小 1.59tVhm 。 年 均 地 表 径流 减少 区 域 占 全 流域 的 89.32% ,主要 位 于 流域 
上 游 的 沽 源 和 围场 ,中游 的 承德 县 和 承德 市 双 桥 区 ,下游 的 青龙 、 建 昌 、 迁 安 与 入 海口 的 三 角 洲 地 区 ;流域 年 均 
产 水 量 增加 区 域 占 全 流域 76.71% , 主要 位 于 流域 上 游 的 正 蓝 族 、 克 什 克 腾 旗 、 多 伦 ,中 游 的 承德 市 双 桥 区 和 宽 
城 ,下 游 的 迁 西 . 青龙 . 迁 安 卢 龙 与 漆 县 ;年 均 产 沙 减 少 区 域 占 全 流域 的 93.89% ,位 于 上 游 沾 源 和 围场 .中游 
承德 县 与 下 游 的 青龙 和 迁 西 。 

(3) 地 表 径流 、 产 水 、 产 沙 与 林地 面积 呈 负 相关 ; 产 水 、 产 沙 与 草地 面积 分 别 旦 正 、 人 负 相 关 ; 地 表 径 流 、 产 沙 
与 农业 用 地 呈正 相关 ;地 表 径流 、 产 水 、 产 沙 与 景观 形状 复杂 程度 、Shannon's HAE ,4) 5 EE 6: YE HS , 53 EE SE 
E .最 大 斑 块 指数 呈 负 相关 ; 产 水 产 沙 与 景观 斑 块 密度 Shannon's 多 样 性 指数 呈正 相关 。 

(4) 京 津 风沙 源 治理 工程 开展 以 来 ,年 均 径流 量 和 产 沙 量 显 著 减少 ,汛期 径流 与 泥 沙 明显 降低 ; 产 水 与 产 
沙 高 值 区 显著 缩小 。 中 上 游 大 部 分 区 域 年 均 产 沙 小 于 1 thm ,大 于 10 t/hm 的 区 域 已 消失 。 产 沙 关 键 源 区 
仍 需 重点 治理 ,应 继续 加 强 林 草 建设 ,提高 植被 质量 与 履 盖 率 ; 对 陡坡 耕地 实施 退耕 还 林 , 对 严重 沙化 耕地 实 
施 退耕 还 草 ;集约 化 用 地 等 ;更 应 重视 有 利于 防治 泥 沙 输出 的 景观 类 型 (如 林地 草地) 的 空间 聚集 程度 、 延 展 
性 .连接 度 与 规则 形状 。 
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